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Abstract 

The BABAR detector, situated at the SLAC PEP-II asymmetric e + e~ collider, has been recording 
data at energies on and around the T(4S) resonance since May 1999. In this paper, we briefly 
describe the PEP-II B Factory and the BABAR detector. The performance presently achieved 
by the experiment in the areas of tracking, vertexing, calorimetry and particle identification is 
reviewed. Analysis concepts that are used in the various papers submitted to this conference are 
also discussed. 
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1 Introduction 



The principal goal of an asymmetric B Factory running at energies on and around the T(4S) 
resonance is the comprehensive study of CP violation in B meson systems. The asymmetry in 
beam energies makes possible the measurement of time-dependent CP violating asymmetries in the 
decay of neutral B mesons. 

CP asymmetries in the B meson system are expected to be large, so that relatively small samples 
of events are needed to provide accurate measurements. Unfortunately, the B meson decay channels 
of interest for the study of CP violation have extremely small branching fractions, of order of 1CP 4 
or below. In order to measure asymmetries with errors at the 10% level or better, samples of several 
tens of million neutral B meson pairs are needed. The new generation of B factories must therefore 
have unprecedented luminosities, in the range 10 33 to 10 34 cm~ 2 s _1 . 

2 The PEP-II B Factory 

The PEP-II B Factory is an e + e~ colliding beam storage ring complex designed to produce a 
luminosity of at least 3 x 10 33 cm _2 s _1 at Ecm = 10.58 GeV, the mass of the T(AS) resonance. 
PEP-II has been constructed by a collaboration of SLAC, LBNL and LLNL on the SLAC site. 
The machine is asymmetric with energies of 9.0 GeV for the electron beam and 3.1 GeV for the 
positron beam. The unequal beam energies require a two-ring configuration: electrons in a High 
Energy Ring (HER) colliding with positrons in the Low Energy Ring (LER). Some accelerator 
parameters and achievements of the High Energy Ring are listed in Table |l|, those of the Low 
Energy Ring in Table The electrons and positrons are produced by the SLAC Linac, whose high 
intensity makes it optimal to recharge the beams in "top-off" mode whenever the luminosity drops 
to about 90% of its peak value. 



Table 1: PEP-II High Energy Ring Performance Results. 



Parameter 


units 


Design 


Best Achieved 


Typical in-run 


Energy 


GeV 


9.0 


9.0, ramp to 9.1 & back 


9.0, ramp 8.84-9.04 


Number of Bunches 




1658 


1658 


553-829 


Total Beam Current 


A 


0.75 (1.0) 


0.92 


0.70 


Beam Lifetime 




4hrs @ 1A 


llhrs @ 0.9A 


9hrs @ 0.70 A 



Table 2: PEP-II Low Energy Ring Performance Results. 



Parameter 


units 


Design 


Best Achieved 


Typical in-run 


Energy 


GeV 


3.1 


3.1 


3.1 


Number of Bunches 




1658 


1658 


553-829 


Total Beam Current 


A 


2.14 


1.72 


1.10 


Beam Lifetime 




4hrs @ 2A 


3.3hrs @ 1.4A 


3hrs @ 1.1 A 



The rings are housed in the 2.2 km former PEP tunnel but with distinct vacuum and accelerating 
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structures. The High Energy Ring (HER) reuses the magnets of the old PEP machine whereas 
the Low Energy Ring (LER) is new and is put in place on top of the HER. The PEP-II design 
has 1658 bunches, each containing 2.1 x 10 10 electrons (HER) and 5.9 x 10 10 positrons (LER), 
spaced at 4.2 ns. The RF system provides a total power of 5.1 MW from seven klystron stations 
driving 24 conventional copper 476 MHz RF cavities. Bunches are brought into a common vacuum 
chamber and septum at a few nTorr and into head-on collisions in an interaction region within a 
2.5 cm-radius beryllium beam-pipe, around which BABAR is located. 

Asymmetric collisions produce a center of mass energy of 10.58 GeV with a boost of /?7 = 0.56 
in the lab frame, thereby allowing a measurement of time dependent CP violating asymmetries. 
Such a boost produces an average separation of fi^cr = 250 fim between the two B vertices, which 
is crucial for studying the cleanest and most promising CP violating modes. 

The PEP-II e + e~ collider became operational in July 1998 with the completion of the LER. 
The first collisions were seen shortly thereafter. Fall and winter 1998 PEP-II runs concentrated on 
raising the beam currents and increasing the luminosity. In February 1999, the peak luminosity 
reached 5.2 x 10 32 cm" 2 s _1 . In a two month spring down time, the BABAR detector was installed. 
PEP-II turned on May W h and BABAR saw its first hadronic event on May 26 th 1999. In August 
1999, PEP-II passed the world record for luminosity, achieving 8.1 x 10 32 cm _2 s _1 . The present 
luminosity in PEP-II is 2.17 x 10 33 cm~ 2 s~ 1 , or about 72% of the design goal. In June 2000, PEP-II 
delivered an integrated luminosity of 150 pb -1 in one day, above the design goal for daily integrated 
luminosity of 135pb~ 1 . Through July 6, PEP-II has delivered over 13.8 fb" 1 to BABAR, of which 
12.7 fb _1 have been logged by the detector. Beam parameters and luminosity performance results 
are shown in Table ||. The present plan is to collide until the end of October 2000 followed by a 
three month installation period. A detailed status report of the current PEP-II performance can 
be found in Ref. |§. 



Table 3: PEP-II Luminosity Performance. 



Parameter 


units 


Design 


Best Achieved 


Typical in- run 


Peak Luminosity 


cm~ 2 s~ x 


3 x 10 33 


2.2 x 10 33 


2.0 x 10 33 


Specific Luminosity 


cm _2 s _1 /bunch 


3.1 x 10 30 


2.9 x 10 30 


1.8 x 10 30 


IP Hor. Spot Size Ux 


filR 


220 


190 




IP Ver. Spot Size Ey 




6.7 


6.0 




Top-off Time 


minutes 


3 


2 


3 


Fill Time 


minutes 


6 


8 


10 


Integrated Luminosity 


pb _1 /8 hours 




54 






pb _1 /day 


135 


154 






fb _1 /week 




0.8 






fb -1 /month 




2.7 




Total Integrated Luminosity 


fb- 1 




15.4 





(thru 7/20/2000) 
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Single Beam Background - HER 




-100 -50 50 100 

Longitudinal Position in Beam Direction (cm) 




10 20 30 40 50 60 70 80 90 100 
day of year 2000 



Figure 1: Left: z distribution of LI pass-through tracks at their point of closest approach to the 
beam line. Right: Time-history of the DCH current, normalized to that expected, at the same beam 
currents and luminosity, from background measurements performed in July 2000. The exponential 
decrease during the first two months of operation reflects the scrubbing, by synchrotron-radiation 
photons, of residual gas molecules off the vacuum-pipe wall. 

2.1 Machine Backgrounds 

Operationally, the acceptable level of backgrounds is determined primarily by the radiation hard- 
ness of the Silicon Vertex Tracker (SVT) and Electromagnetic Calorimeter (EMC) detectors, and 
by requiring a tolerable Drift Chamber (DCH) current. The Level-1 (LI) trigger rate and the 
occupancy in the other detector systems also constitute occasional limitations. Careful measure- 
ment, analysis and simulation of the background sources and their impact have led to a detailed 
understanding and an effective remediation of their effects. 

The primary causes of steady-state backgrounds in PEP-II are, in order of increasing impor- 
tance: 

• Synchrotron radiation generated in the bending magnets and final focusing quadrupoles in 
the incoming HER and LER beam lines. Careful layout of the interaction-region area and a 
conservative synchrotron radiation masking scheme have proven very effective against these 
sources. 

• The interaction of beam particles with residual gas around the rings (beam-gas) , which con- 
stitutes the primary source of radiation damage and has had, averaged over this first run, the 
largest impact on operational efficiency. 

• Collision-related electromagnetic shower debris, dominated by energy- degraded from 
radiative-Bhabha scattering which strike vacuum components within a few meters of the 
interaction point (IP). This background, directly proportional to the instantaneous luminos- 
ity, was barely detectable in early running; it now noticeably affects all detectors except the 
SVT. 
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While instantaneous background conditions vary because of interaction point orbit drifts and 
sensitivity of beam tails to small tune adjustments, reproducible patterns have emerged. 

The relative importance of the single-beam and luminosity contributions varies from one detec- 
tor system to the next. The HER beam-gas contribution significantly impacts all detectors except 
the particle identification system (DIRC): the combination of a 40 m long straight section, almost 
devoid of magnetic bending upstream of the final doublet, and of the magnetic beam- separation 
scheme, result in abundant bremsstrahlung-induced electromagnetic debris being directed onto the 
IP vacuum pipe. The same process occurs in the LER, but to a lesser extent (for a given beam cur- 
rent) because of a shorter drift section and lower primary energy. Most detectors, therefore, exhibit 
occupancy peaks at (j> = and 180°, reflecting the fact that the separation dipoles bend energy- 
degraded particles in opposite directions. Such local beam-gas interactions contribute noticeably 
to the DCH current; they dominate the SVT instantaneous dose rate in the horizontal plane, the 
EMC crystal occupancy above the cluster-seeding threshold (10 MeV), and the LI trigger rate. 
Analysis of pass-through events from the charged LI trigger shows that the beam-pipe wall and 
several aperture restrictions within 100 cm of the interaction point are the primary impact points 
for lost particles, as illustrated for the HER in Figure |l| (left)). To minimize this background it is 
vital to maintain a low pressure in the region from 4 m to 60 m upstream of the interaction point 
in the incoming HER and LER beam lines. Scrubbing, which has reduced the average dynamic 
pressure in both rings to levels close to or below the design value, also proved effective (Figure |l| 
(right)). 

Under typical running conditions, the DIRC and the DCH receive comparable contributions 
from beam-gas Coulomb scatters around the entire LER, and from luminosity backgrounds. In 
addition, these two detectors proved particularly vulnerable to tails generated by beam-beam or 
electron-cloud induced blowup of the low-energy beam. Even though partially eliminated by a set 
of betatron collimators in the last arc, such tails - which are very sensitive to accelerator conditions 
and tuning procedures - tend to scrape near the highest-/3 point of the final LER doublet, located 
inside the DIRC's standoff-box (SOB). The resulting background fluctuations occasionally cause 
the DCH high-voltage protection system to trip, or the counting rate in the DIRC photomultipliers 
to exceed 200 kHz. The problem could be only partially alleviated by local lead shielding; additional 
collimation will be installed during the fall 2000 shutdown. 

While trigger-rate and occupancy considerations define acceptable dynamic running conditions, 
it is the total integrated radiation dose that determines the lifetime of the detector systems. Despite 
a significant investment in radiation-hard technology, the innermost silicon layers of the SVT and 
its front-end electronics remain the most susceptible to radiation damage. The accumulated dose 
has been maintained below budget, as shown in Figure [2|, through a strict program of hardware 
interlocks, administrative controls, and real-time monitoring. As of this conference, the worst 
irradiated spot of the SVT has been exposed to approximately 300 kRad, 20 to 30% of which is 
contributed by injection periods. 

3 The BABAR Detector 

The new magnetic spectrometer BABAR (Figure ||) has been constructed at SLAC, by a collaboration 
of nine countries, to precisely measure e + e~ annihilations at center of mass energy y/s ~10GeV 
produced with the new PEP-II Asymmetric Storage Rings. Construction of this detector was 
approved in November 1995 and its commissioning followed in the Fall of 1998. First data with 
PEP-II colliding beams was collected in May 1999. 
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Figure 2: projected History plot of accumulated radiation dose in the SVT using PIN diodes near 
the interaction point in the horizontal plane (the black line represents the allotted budget, rationing 
out a 2 MRad dose over the detector's planned life span). 

The BABAR superconducting solenoid, which produces a 1.5T axial magnetic field, contains a 
set of nested detectors: a five-layer silicon-strip vertex detector (SVT), a central Drift Chamber 
(DCH) for charged particle detection and momentum measurement, a quartz-bar Cherenkov radi- 
ation detector (DIRC) for particle identification, and a Csl crystal Electromagnetic Calorimeter 
(EMC) for detecting photons and separating electrons from charged pions. The calorimeter has a 
barrel section, and an endcap which extends it asymmetrically into the forward direction (e~ beam 
direction), where many of the collision products emerge. (There is no calorimeter coverage in the 
backward direction.) Two layers of cylindrical resistive plate chambers (RPCs) are located between 
the barrel calorimeter and the magnet cryostat. All the detectors located inside the magnet have 
full acceptance in azimuth. The Instrumented Flux Return (IFR) outside the cryostat is composed 
of 18 layers of steel, which increase in thickness outward, and contains 19 layers of planar RPCs in 
the barrel and 18 in the endcaps. The RPCs allow the separation of muons and charged hadrons, 
and also detect penetrating neutral hadrons (mainly Kl). 

The following sections provide more details on the subsystems of the detector together with 
current performance results. 

3.1 Silicon Vertex Tracker 

The Silicon Vertex Tracker (SVT) provides the required vertex resolution for the measurement 
of CP violation at PEP-II. In addition, it is capable of independent charged particle tracking, of 
particular importance for low transverse momentum particles (pr < 120MeV/c) that cannot be 
measured by the central tracking chamber. The detector design is optimized to take into account 
the physical constraints imposed by the PEP-II geometry at the interaction region: the presence of 
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Figure 3: The BABAR Detector. 1. Silicon Vertex Tracker (SVT), 2. Drift Chamber (DCH), 
3. Particle Identification Subsystem (DIRC-Detector of Internally Reflected Cherenkov Light, 4. 
Electromagnetic Calorimeter (EMC), 5. Magnet, 6. Instrumented Flux Return (IFR). 

the Bl permanent magnets at ±20 cm from the interaction point, which separate the beams after 
head-on collisions. The acceptance in polar angle 9 is limited by the gap between beamline elements 
to -0.87 < cosOiab < 0.96 (-0.95 < cos6> cm < 0.87). As the innermost BABAR subdetector, two 
important considerations in optimizing the design were low mass, to minimize multiple scattering, 
and radiation hardness of its components. A detailed description of the SVT and its components 
can be found in Ref. Q. 

The SVT contains 52 modules built out of AC-coupled double-sided silicon sensors (strips 
othogonal on the two sides) . These are read out by a full-custom low- noise radiation-hard integrated 
circuit, known as the AtoM chip (mounted on a passive hybrid circuit that is attached to a water 
cooling channel). The detector modules are organized in five radial layers, each divided azimuthally 
into 6, 6, 6, 16 and 18 staves respectively (see Table ^|). The three inner layers are crucial for vertex 
and tracking resolution, while the outer two layers are needed to provide additional measurements 
for stand alone tracking independent of drift chamber information. These outer two layers each 
contain two different types of modules, an inner (labelled a in the Table Q) and an outer (labelled 
b) layer, occupying slightly different radial positions. The modules are assembled on carbon fiber 
support cones, which in turn are positioned around the beam pipe and the Bl magnets. The SVT 
and some beamline elements are housed inside a strong support tube, with its load transferred at 
the ends to the PEP-II beamline support "rafts." 
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SVT Hit Resolution vs. Incident Track Angle 
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Figure 4: SVT single point resolution as a function of the incident track angle for the innermost 
layer. 



Table 4: Layer Structure of the SVT. 
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A system consisting of 12 PIN photodiode sensors is placed close to the first silicon layer to 
monitor continuously the radiation exposure of the SVT and manage its radiation dose budget. 

During the first year of data taking, all major design goals for the detector were already achieved. 
The average hit reconstruction efficiency is above 98% in both views. The hit resolution presently 
observed, following an initial accurate alignment procedure, is shown in Figure ||. This compares 
well with Monte Carlo predictions for the resolution with perfect alignment. The impact parameter 
resolution, dominated by the SVT's precision for measuring high transverse momentum tracks, is 
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shown in Figure ||. Two-track vertices, such as J/ip — > /^ + ^ _ are reconstructed with a typical 
resolution of 50 /im; the resolution on the z separation between the two B decay vertices is typically 
110 /xm, in good agreement with the design goals. 
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Figure 5: Impact parameter resolution in the r — <f> (ado m the left plot) and in the z-view (a z in 
the right plot) as measured on data. The horizontal line in the left-hand plot is 1/V2 times the 
transverse miss distance for muon pair events. 



3.2 Drift Chamber 

The Drift Chamber (DCH) is the main tracking device for charged particles with transverse mo- 
menta pt above about 120 MeV/ c, providing a precision measurement of px from the curvature of 
charged particle trajectories in the 1.5 T magnetic field. Prompt single-cell hit information from 
the chamber is also a basic component of the Level-1 trigger (see Section 3.6.1). The DCH is a 
280 cm-long cylinder, with a inner radius of 23.6 cm and and outer radius of 80.9 cm. It is mounted 
in cantilever from the rear endplate to the DIRC central support tube, within the volume inside 
the DIRC and outside the PEP-II support tube. The center of the chamber is displaced forward by 
36.7 cm to improve the forward track length, given the asymmetric boost for the T(4£) events. The 
beryllium inner wall (0.28% radiation lengths), the thin outer half of the forward endplate (15 cm 
aluminum), and the carbon- fiber outer cylinder are all designed to minimize material in front of 
the calorimeter. 

The drift system consists of 7104 hexagonal cells, approximately 1.8 cm wide by 1.2 cm high, 
arranged in 40 concentric layers between a radius of 25.3 and 79.0 cm. The active volume provides 
charged particle tracking over the polar angle range —0.92 < cos 9i a b < 0.96. The forty layers are 
grouped into ten superlayers of four layers each, organized with the same orientation for sense and 
field wires within a given superlayer. The superlayer structure facilitates fast local segment finding 
as the first step in pattern recognition. This arrangement is particularly important for Level-1 
trigger decisions. Superlayers alternate in orientation: first axial (A), then a small positive stereo 
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angle (U), followed by a small negative stereo angle (V). All superlayers participate in the Level-1 
trigger track finding; only the axial superlayers participate in the Level-1 trigger pt determination. 



Table 5: Layer arrangement for the DCH 



Superlayer 


Inner 


Cells/ 


Stereo 




Radius [cm] 


Layer 


Angle [mr] 
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26.04 


96 
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31.85 


112 


+ [44.9-50.0] 


3 


37.05 


128 


-[52.3-57.4] 


4 


42.27 


144 





5 


48.08 


176 


+ [55.6-59.7] 


6 


53.32 


192 


-[62.8-66.9] 


7 


58.54 


208 





8 


64.30 


224 


+ [65.0-68.5] 


9 


69.52 


240 


-[72.1-75.8] 


10 


74.72 


256 






The individual hexagonal cells each consist of a 20 /mi rhenium-tungsten sense wire operating 
nominally in the range 1900-1960 V, surrounded by 6 cathode wires, approximately half of which 
are shared with neighboring cells. For the inner cells in a superlayer, the cathode wires are 120 /im 
aluminum, grounded to the rear endplate. For outer cells one of the cathode wires is 80 fim 
aluminum, held at about 350 V. All sense and cathode wires are gold plated. Table |H| summarizes 
the main geometric features of the design. The counting gas, chosen to have low mass, is an 
80%: 20% mixture of Hedsobutane, with a small amount (3000 ppm) of water vapor to prolong 
the lifetime of the chamber in a high rate environment. The combination of low- mass gas and 
minimization of material in the field cage itself is designed to reduce the contribution of multiple 
scattering to the px resolution for the typically soft tracks produced in T(4S) events. 

The DCH readout electronics, including front-end amplifier, timing and pulse height digitiza- 
tion, event readout and high voltage distribution, are all mounted on the rear endplate (opposite the 
boost direction) to keep the additional material outside the detector fiducial volume. The instru- 
mentation includes two custom integrated circuits: a low-noise bipolar amplifier IC and a CMOS 
digitizer IC that incorporates 8 channels of TDC, Flash-ADC, pipelined data-storage, multiple 
event buffers and prompt trigger output. Data readout is multiplexed on four fiber-optic links to 
BAR4R-standard VME readout modules (ROM); trigger data is transported on 24 fiber-optic links. 
The ROMs extract the integrated charge from the digitized waveform before passing the data on to 
next level in the data acquisition system. High voltage is supplied by a CAEN SY527 mainframe 
and is distributed through circuit boards that plug directly on the DCH endplate feedthroughs. 
Temperature and humidity sensors monitor the chamber environment, and radiation monitors track 
the accumulated dose. 

The nominal design for the chamber calls for an average 140 /jm single point resolution. The 
space-time relation (STR) for the non-saturated counting gas has been modeled with a separate 
6* ft -order Chebychev polynomial for the left and right parts of the cell. The STR has proven to be 
very stable over time. There is a small residual dependence on gas density, not yet corrected in 
production. The single cell resolution function obtained from the ensemble of charged tracks in the 
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normal event stream is shown in Figure |6] for operation at 1960 V. This shows a weighted average 
resolution of 125 /im. 

In addition to recording timing information for hits, the Flash- ADC is used to measure the time 
development of pulse height in the cell. Signal processing (feature extraction) in the ROM converts 
this into a measurement of deposited energy. Real-time overall gain corrections are calculated and 
applied during the reconstruction step. Saturation, path length, polar angle, individual wire and 
layer gain corrections are applied as well, although these corrections are stable over long periods 
of time. The resulting distribution of the truncated mean dE/dx observed for Bhabha electrons is 
shown in Figure |6| The observed resolution is found to be 7.5%. We expect to achieve the predicted 
resolution of 7% after additional corrections, including the effect of varying entrance angles. 
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Figure 6: DCH single cell resolution (left) and dE/dx resolution for Bhabha electrons (right). 



3.3 Particle Identification System 

Standalone identification of charged particles is based on a specialized subdetector system that uses 
the Detection of Internally Reflected Cherenkov light (DIRC). Charged particles exiting the barrel 
region of the DCH transit an array of 144 fused silica quartz bars, each approximately 17 mm thick 
(5r), 35mm wide (5(r<p)) and 4.9m long (Sz). Particles above Cherenkov threshold radiate photons 
in the quartz media. The angles of the Cherenkov photons with respect to particle direction are 
measured with an array of 10,752 photomultiplier tubes located outside the return yoke, in a special 
low magnetic field volume. The polar angle coverage is —0.84 < cos 9i a b < 0.90. 

The 144 quartz bars of the DIRC are arranged in 12 modules, or "bar-boxes", that penetrate 
through the magnetic end plug. Cherenkov photons trapped in a quartz bar exit the bar through 
a wedge and a quartz window into a water tank, which optically couples the quartz bars to the 
photomultiplier array. The photomultiplier tubes are arranged in 12 sectors corresponding to the 
12 barboxes, on the surface of a half torus, with a major radius of about 3 m and a minor radius 
of 1.2 m. By knowing the particle direction from the tracking system, and the location of the 
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photomultiplier tube that observes a Cherenkov photon, the angle that a Cherenkov photon makes 
with respect to the track direction can be deduced. Due to internal reflections within the bars, 
there are several solutions for a hit-to-track association. The angular resolution for a single photon 
is about 10.2 mr (see Figure 0(a)) and, with an average of 30 photons per track, the "per track" 
Cherenkov resolution is about 2.8 mr, rms. This corresponds to a separation of approximately three 
standard deviations between charged ICs and 7r's at 3GeV/c. An additional observable is the time 
of a photomultiplier hit, measured to a precision of 1.7ns (see Figure ^(b)). By comparison of 
the measured hit time with the propagation time of a photon solution, it is possible to effectively 
suppress background photons and invalid photon solutions. 




Figure 7: Resolution of the reconstructed Cherenkov angle for single photons (a) and of the differ- 
ence between measured and expected arrival time (b). 



3.4 Electromagnetic Calorimeter 

The Electromagnetic Calorimeter (EMC) contains 6580 Csl crystals doped with thallium iodide 
at about 1000 ppm. Each crystal is a truncated trapezoidal pyramid and ranges from 16 to 17.5 
radiation lengths in thickness. The front faces are typically ~5cm in each dimension. The crystals 
are arranged quasi-projectively in a barrel structure of 48 polar {0) rows by 120 crystals in azimuth 
((/)), with an inner radius of 90 cm. The forward end is closed by a separable endcap capable of 
holding nine additional rows. At present, eight are filled with crystals and the innermost with 
shielding absorber. The polar angle coverage of the calorimeter is —0.78 < cos 9^ < 0.96. Beam- 
line elements occlude cos 9i a b > 0.94. Each crystal is wrapped with a diffuse reflective material 
(TYVEK) and housed in a thin eggcrate-like carbon fiber composite mechanical structure. There 
are 280 such modules in the barrel (7 types, 40 of each type) and 20 identical endcap modules. 
Crystals are read out with two independent 2 cm 2 large area PIN photodiodes epoxied to their rear 
faces. Dual-range ASIC preamplifiers reside directly behind the photodiodes in a shielded housing 
that also provides a thermal path for heat removal. Shielded ribbon cables carry analog signals 
to the end flanges of the barrel and the back plate of the endcap, where additional amplification 
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and digitizing electronics provide for a total of four overlapping linear ranges. The system handles 
signals from ~50keV to ~13GeV, corresponding to 18 bit dynamic range. A short shaping time 
of ~ 400 ns is used in the preamplifiers to reduce the impact of soft (<5MeV) beam-related pho- 
ton backgrounds. Noise performance can be recovered by digitally processing the signal waveform 
sampled at 4 MHz. Calibration and monitoring is achieved by charge injection into the front end 
of the preamplifiers, a fiberoptic/xenon pulser system injecting light into the rear of each crystal, 
and a circulating radioactive source (a neutron-activated fluorocarbon fluid) producing a 6.13 MeV 
photon peak in each crystal. Signals from data (vr°s, rjs, radiative, and non-radiative Bhabhas, 
77 and fj,^ yT events) can provide additional calibration points. Source and Bhabha calibrations 
are updated weekly to track the small changes in light yield with integrated radiation dose. Light 
pulser runs are carried out daily to monitor relative changes at the <0.3% level. 

The calorimeter has achieved an electronics noise energy (ENE) of ~220keV (coherent plus 
incoherent) measured with the source system (in the absence of colliding beams) after digital signal 
processing. During regular data taking, this digital filtering is not applied and the ENE rises to 
~450keV owing to the short shaping time; consequently, only channels with >1 MeV are presently 
used in the reconstruction of calorimeter energy deposits. The efficiency of the calorimeter exceeds 
96% for the detection of photons with energies above 20 MeV. 

The energy resolution can be measured directly with the radioactive source at low energies and 
with electrons from Bhabha scattering at high energies, yielding resolutions of <r(E)/E = 5.0 ±0.8% 
at 6.13 MeV and <r(E)/E = 1.9 ± 0.07% at 7.5 GeV, respectively. The energy resolution can also 
be inferred from the observed mass resolutions for the tt° and rj, which are measured to be around 
7 MeV and 16 MeV, respectively. 
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Figure 8: The energy resolution as a function of energy, as determined from the observed width 
of 7T° and 77 decays to two photons of equal energy, and the resolution for Bhabha electrons. The 
shaded band is the best fit to the 7r°, r], and Bhabha data. Also shown is the energy resolution of 
the 6.13 MeV photons from the radioactive source, and of the photons in the transition \ci ~~ * J/^l- 
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Figure || shows the energy resolution extracted from a variety of data as a function of energy. 
A fit to the 7T°, r], and Bhabha energy resolution measurements, assuming an energy dependence 
of the form: 

<t(E)/E = ai(E/GeVr 1/4 ©a 2 , (1) 

gives <7i = (2.32±0.30)% and a 2 = (1.85±0.12)%. Implementation of digital filtering will reduce the 
electronics and beam background noise. This, together with the afforded improvements in cluster 
reconstruction algorithms, will lead to better energy resolution for the 2001 run, particularly at 
lower energies. 

The 7T° and 77 data are also used to measure the angular resolution of the calorimeter. It is 
found to vary between about 12 mr at low energies and 3mr at high energies, described by an 
energy dependence according to: 

ae, (t> = a 1 (E/GeV)- l/2 + a 2 , (2) 

with at = (3.87 ± 0.07) mrad and a 2 = (0.00 ± 0.04) mrad. 

The calorimeter is also used for the separation of hadrons from electrons, and in conjunction 
with the IFR for muon and K® identification; its performance in these areas is presented in section 
7.2. 

3.5 Instrumented Flux Return 

The Instrumented Flux Return (IFR) is used to identify muons and neutral hadrons. The detector 
consists of nearly 900 Resistive Plate Chambers (RPCs) interleaved with the iron plates which 
comprise the flux return for the 1.5 T solenoidal magnet. There are 19 RPC layers in the barrel 
region and 18 layers in the forward and backward endcaps. The iron plate thickness is graded from 
2 cm for plates closest to the interaction region to 10 cm for the outermost layers for a total of 
>65cm in the barrel and >60cm in the endcaps. An additional cylindrical RPC containing two 
layers is located in the barrel region between the calorimeter and solenoid cryostat. 

The single gap planar RPCs || are constructed from 2 sheets of bakelite separated by 2 mm 
spacers. The inner surfaces are treated with linseed oil. The outer bakelite surfaces are painted 
with graphite. Layers of insulating mylar, aluminized mylar pickup strips, foam and an aluminized 
mylar ground-plane are glued to both sides of the bakelite forming the RPC package. The pickup 
strips are connected with discriminator cards which record a bit per strip if an in-time signal above 
threshold is recorded. The pickup strips on top and bottom are orthogonal, providing x-y position 
information in the endcaps and 4>-z information in the barrel. The cylindrical RPCs are constructed 
from a special plastic requiring no surface treatment. The RPCs utilize a gas mixture of 56.7% 
Argon, 38.8% Freon-134a, and 4.5% Isobutane and operate at 7.6 kV. The cylindrical RPCs use the 
same gas but run at 7.4 kV. The steel flux return is water-cooled to reduce temperature excursions, 
keeping the barrel at ~20°C and the endcaps at ~22-24°C 

The RPC efficiencies are measured in cosmic ray runs taken weekly and in collision data. The 
average chamber efficiencies during the 2000 run were ~78% in the barrel, and ~87% in the forward 
endcap, lower than the average ~92% measured in June 1999. A small fraction of the RPCs is 
presently disconnected because of high-voltage problems. A vigorous testing and R&D program is 
underway to understand the efficiency decrease and to improve RPC performance. 
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3.6 Trigger 

The BABAR trigger system consists of the Level- 1 hardware trigger and the Level-3 software trigger. 
The Level-1 trigger decision is issued within a latency of 11-12 ^s after the corresponding beam 
crossing to initiate the Data Acquisition (DAQ) readout of the relevant time slice of the pipelined 
data for all detectors. The latency for hadronic events is confined to a smaller, ±150ns, range. 
The acquired data are processed by Level-3 running on the online farm nodes to select physics and 
calibration events for data logging. 

3.6.1 Level-1 Trigger 

The Level-1 (LI) trigger system comprises the DCH track trigger (DCT), the EMC energy trigger 
(EMT), the IFR muon trigger (IFT) and the global trigger (GLT). The LI trigger system operates 
in a continuous sampling mode, processing input data and generating output trigger information 
at fixed time intervals. The main logic of the LI trigger is implemented in six types of 9u VME 
modules with a total of 41 boards. The basic trigger primitives generated by the DCT, EMT and 
IFT, mostly cf> maps of the particle signatures, are collected by the GLT to form 24 trigger lines for 
the fast control system. The DCT identifies short and long track primitives as well as high pt track 
primitives. Using the loosest criteria, short tracks are accepted down to pr > 120MeV/c if they 
reach superlayer five in the DCH. The EMT generates trigger primitives passing various energy 
thresholds. The lowest threshold is 100 MeV, in order to be fully efficient for minimum ionizing 
particles. 

All the LI trigger system components appear to the data acquisition system as standard detector 
subsystems, in addition to their role in providing the trigger signals. The data read out from them 
provides details on the operation of the trigger for each acquired event, facilitating diagnostics and 
the determination of the LI trigger efficiency. 

The LI trigger system is designed to be able to trigger independently from pure DCT and pure 
EMT triggers with high efficiency for most physics sources, in order to allow cross-calibration of 
efficiencies. In particular, BB events are triggered at > 99% efficiency from either DCT or EMT 
and the combined efficiency is > 99.9%. Only r and two-photon events do not have a fully efficient 
pure EMT trigger and rely mainly on DCT triggers. 

The total LI trigger rate for a typical run with HER (LER) currents at 700 mA (1100 mA) and 
a luminosity of 2. Ox 10 33 cm~ 2 s _1 is ~700 Hz. This is typically stable to within 10% for the same 
machine currents but can be 50% higher after a major shutdown when the vacuum is relatively 
poor. The joint trigger/data acquisition system rate capability is comfortably beyond the design 
rate of 2 kHz. This avoids significant deadtime, and leaves room for the expected increase in 
luminosity. Within the typical LI trigger rate of 700 Hz at 2xl0 33 cm _2 s _1 , Bhabha events and 
other e + e~ interactions amount to around 120 Hz, and cosmic ray interactions to 130 Hz. The 
dominant source of background causing the remaining triggers is the interaction of lost particles 
with the beam line components. The LI background rate coming from the HER beam is ~3 times 
higher than that from the LER beam for the typical operating currents. 

3.6.2 Level-3 Trigger 

The Level-3 (L3) trigger is the first component of the trigger system to see complete events. It is 
embedded in the Online Event Processing (OEP) framework running in parallel on 32 event filter 
farm nodes. The L3 trigger processes data from the DCH track trigger and from the EMC using 
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two independent algorithms that form track and cluster objects. These are fed into a set of filter 
algorithms that operate on either one or both. The L3 DCH algorithm performs fast lookup-table 
based track finding and three-dimensional track fitting, efficient for tracks with pr > 250 MeV/ c. 
The event to is determined from DCH trigger track segments and hit data to better than 5 ns. The 
L3 EMC algorithm uses calorimeter crystal data in a fast two-dimensional clustering algorithm 
implemented using a lookup table. To reduce noise, the crystal energy threshold is set to 30 MeV. 
L3 clusters are accepted for energies > 100 MeV. 

The L3 logging decision is based on generic track/cluster topologies rather than on the iden- 
tification of individual physics processes. An exception is made for Bhabha events, which have to 
be vetoed to reduce their rate (~ 100 Hz at 2 x 10 33 cm~ 2 s~ 1 luminosity). The physics trigger is a 
logical OR of two orthogonal filters. The track filter requires either one track with pr > 800MeV/c 
coming from the interaction point, or two tracks with px > 250MeV/c and slightly looser vertex 
cuts. The cluster filter accepts events with a high multiplicity or with a large total energy in the 
EMC and a high invariant (pseudo) mass. The efficiency of the track filter for BB events is 99%, 
while the efficiency of the cluster filter is 94%. 

Both filters are subject to a veto algorithm that identifies Bhabha events based on clean signa- 
tures in the DCH and the EMC. The veto has no impact on hadronic events, and only has noticeable 
effect on very few types of events such as rr — > evvevv. A prescaled sample of Bhabha events, 
flattened in 6, is preserved for calibration purposes. In addition, L3 logs various other samples 
for calibration and monitoring, such as radiative Bhabha events, cosmics, and random triggers. 
The luminosity is measured online by L3 using a track-based algorithm with a scale calibration 
precision of <~ 5%. The typical logging rate at 2.0 x 10 33 cm _2 s _1 is ~ 85 Hz, well within the 
design specifications. 

3.7 Data Acquisition and Online System 

The BABAR data acquisition system employs a set of standard Readout Modules (ROMs) to commu- 
nicate with the front-end electronics of the detector systems. A common communications fabric — 
1 Gbps optical fiber links — and protocol is used for readout by all the systems, with detector-specific 
extensions to the protocol for control. 

Each ROM is built around a commercial VME single-board computer running Vx Works, linked 
with three custom boards, one of which supplies one of two "personalities". The calorimeter uses 
an untriggered version of the ROM which continuously collects data from up to three input data 
fibers, separating out in software the data corresponding to LI triggers. All of the other detector 
systems use a version which requests and collects event data from up to two fibers upon the receipt 
of a LI trigger. 

ROMs are organized into crates in groups of two to eleven, including a master ROM. The master 
builds partial events from the data acquired by the crate's other ROMs from the front ends. Final 
event building is performed by sending these data fragments via 100-Mbps switched Ethernet to 
the nodes of the Online Event Processing (OEP) farm, which are commercial Unix workstations. 
Static load-balancing is used, allocating events to nodes using a simple deterministic algorithm 
based on a unique time stamp assigned to each LI trigger. 

The farm nodes in OEP apply the L3 trigger algorithms and perform fast data quality monitor- 
ing functions. Events passing the L3 selection are sent via TCP/IP to a single server process which 
logs them to disk files. Data quality monitoring results are collected from all nodes and merged for 
display to operators and for automated comparison against defined references. 
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This system has easily exceeded its fundamental performance requirement of acquiring 2000 Hz 
of LI triggers and logging L3 triggers at up to 100 Hz. In recent running we have been operating 
far from its limits, with a LI rate of approximately 700 Hz, an average event size of 27 kB, and a 
L3 logging rate of 85 Hz. There is therefore considerable headroom already available for operation 
at higher luminosities, but the system can be scaled further, if necessary, up to the data bandwidth 
limits of the front-end electronics through the addition of more ROMs and more OEP nodes. 

The online computing system includes several additional components. The Detector Control 
system, which monitors and controls the experiments environmental systems such as low and high 
voltage power and detector gas supplies, as well as monitoring parameters of the PEP-II collider. 
The Run Control subsystem coordinates the activities of all the online and data acquisition com- 
ponents, and interlocks data acquisition with the maintenance of safe and acceptable detector and 
collider conditions. 

3.8 Online Prompt Reconstruction 

The first stage of production processing after data acquisition, the Online Prompt Reconstruction 
or OPR, is performed automatically by a compute farm which currently consists of about 150 Unix 
processors. All colliding beam events are filtered and tagged. Interesting events are completely 
reconstructed and written to the object database. A novel system of "rolling calibrations" extracts 
updated reconstruction constants from the data which are then written to the object database and 
used for subsequent reconstruction. Finally, detailed monitoring distributions are generated and 
fed back into the experiments data quality monitoring program. 

OPR is designed and sized to keep up with the data acquisition with minimum latency. Typ- 
ically, this latency between data availability and processed data is in the range from one to four 
8-hour shifts. Thus, the system operates 24 hours per day and seven days per week. 

A second, identical farm of processors is used for data reprocessing. Using the same mechanisms 
as OPR, the reprocessing system is used to reprocess older data once newer code or improved 
constants become available. 

3.9 Computing System 

The BABAR computing system includes the online system components described above, plus the 
computing required for offline reconstruction, simulation and analysis. This is implemented on 
several hundred commercial Unix computers (currently running Solaris, Linux and OSF operating 
systems), functionally distributed, with tcp/ip-based networking. Fileservers with RAID disk arrays 
and access to HPSS-based mass tape storage in STK robots are used to store and serve the data. 
Network switching employs Cisco technology both online and in the SLAC computing center. 

The BABAR software is mostly C++, designed with OOAD methodology. A commercial database, 
Objectivity, is used for data and condition storage. Code distribution employs BABAR-developed 
release tools built upon CVS and the AFS distributed file system. The LSF batch scheduling 
system is used at the SLAC site. Data distribution to remote sites uses a combination of WAN, 
DLT7000 and Redwood/Eagle tapes. Data acquisition and reconstruction is performed at SLAC. 
Simulation and data analysis tasks are shared between SLAC and remote sites. These remote sites 
include major centers in France, Italy, UK and USA. 
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4 Data Sample 



Results for this Conference are based on data collected up to July 6, 2000, when the total integrated 
luminosity recorded by the experiment was 12.71b" 1 . Ofthis sample, 11.31b" 1 was taken "on peak" , 
near the peak of the T(4S) resonance, and 1.4 fb -1 was taken "off peak" at a center-of-mass energy 
40 MeV below the resonance. For the Osaka conference, the physics analyses are based on a sample 
of 8.9 ftT 1 , of which 7.7 fb -1 is on peak and 1.2 fb 1 off peak. An additional 1.1 fb 1 of on-peak 
data was used in the analysis for sin2/3, for a total of 10.0 fb" 1 . 

The remaining 2.7 fb -1 of recorded data falls into two categories. Data taken during the period 
between June and October 1999 (1.31b" 1 ), when the DIRC system was incomplete (only 5 of the 
12 bar boxes were installed) and various subsystems were still optimizing their operation, will be of 
limited use for future physics analyses. Some data (1.2 fb" 1 ) was excluded from the Osaka sample 
because optimal alignment constants were unavailable when it was processed initially. This data 
will be available after reprocessing in August 2000. 



4.1 Delivered and Recorded Integrated Luminosity 

The luminosity is monitored by PEP-II and by BABAR in the trigger and prompt reconstruction. 
The L3 trigger calculates the deadtime, which is typically less than two percent. It also measures 
the luminosity from Bhabha events. The offline determination (described in Section 8.4) uses 
dimuon and two-photon events as well for cross-checks. The DAQ efficiency is monitored closely 
and is typically greater than 95%. Occasional scans of the resonance are used to make sure that 
the experiment is running on the T(4S) peak. 



4.2 Off-Resonance Running 

By running below the T{4S) resonance, BABAR can study the non- BB contribution on-peak. This 



is used to count the number of observed BB (see Section 8J3) and for background subtraction 
in many analyses, especially inclusive measurements. The amount of luminosity devoted to off- 
peak running is of less importance for clean channels like B° — ► J/t/jK®. The uncertainty of a 



measurement with a ratio of signal plus background to signal b depends [12| on the fraction of 
continuum luminosity c running as: 

~ Tr : ■ (3) 



V <1 - c) 

For an analysis with b = 0.02, collecting about 12% of luminosity off-peak is optimal. This was the 
goal for the year 2000 run, which concentrates on measuring sin2/3. As more luminosity becomes 
available, a somewhat higher rate of continuum running may be chosen to improve the background 
subtraction for modes like B° —* 7r + 7r". 

Operationally, the off-peak running took place about once a month. This minimized the amount 
of time spent tuning the machine after changing the beam energy, while still allowing the analyses to 
track changes in machine backgrounds, reconstruction code, and detector performance. The center- 
of-mass energy was reduced by lowering the HER energy by about 70 MeV, giving a center-of-mass 
energy of 10.540 GeV. The LER energy was left unchanged for simplicity's sake; the resulting 
decrease in boost was negligible for all analyses. 
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4.3 Data Quality 



Data quality is checked at each stage from the trigger to the analysis skims. Gross hardware 
problems, such as power supply failures or trips, are rare and are easily identified by the extensive 
slow control system. A member of the shift crew is assigned to monitor histograms, which are 
updated live with occupancy, timing, event size, and other performance information from the LI and 
L3 triggers, the detector subsystems, and special online event processing cross-checks. Automatic 
processes test additional live histograms and report anomalies to the shifter or subsystem expert. 
This insures that problems are identified and fixed as soon as possible. Several subsystems also 
employ automated processes that perform tests on the data files to make long-term strip charts of 
the kind of information listed above. 

The reconstruction job contains many histograms for each subsystem (including particle iden- 
tification, beamspot and alignment monitoring, and results from some basic analyses). A subset is 
inspected by subsystem experts for each run. The most important quantities are extracted from 
the histograms to make charts showing the long-term behavior of the detector and reconstruction 
code, allowing detection of changes that are too subtle or gradual to notice on a run-by-run basis. 
The quality assessments from live data taking, online processing, and first-pass offline analysis are 
combined to define lists of good runs for further analysis. The yield per unit luminosity versus run 
for each skim is inspected to insure that the data used in physics analyses is of consistent quality 
within samples. Less than one percent of data was declared unusable by this process. 

5 Simulation 

The full detector simulation is structurally divided into three parts: event generation, the tracking 
of particles through the detector, and the detector response simulation. 

5.1 Event Generators 

The simulation process begins by generating events using one of several possible event generators. 
The beam energies used in the simulation are thrown using a Gaussian distribution of width 5.5 
and 2.7 MeV for the high and low energy beams respectively. A transverse momentum spread 
of 2.8 and 1.3MeV/c, obtained from beam emittances for the HER and LER respectively, is also 
simulated. The position of the interaction point is modeled as a Gaussian distribution having 
widths of 125 /im, 4.2 /im, and 0.85 cm, in the x, y, and z directions respectively. The direction 
of the simulated beam axis is rotated with respect to the BABAR detector coordinate system by 
— 18.8mr in the y direction, as observed in data. 

The main generator for BB physics is called EvtGen, details of which can be found in Q. This 
generator provides a framework in which specific decay models can be implemented as modules. 
These modules can perform a variety of different functions such as calculating amplitudes for the 
decay. EvtGen introduces mixing by generating decays of the T(45) to the proper mixture of B°B°, 
B°B°, and B°B° final states, with the correct distributions in At. CP asymmetries are generated 
in modules that modify the generated lifetime distributions of the two .B's produced in the decay of 
the T(4S). Generic models are available for simulating two-body decays to a pair of scalar mesons, 
a scalar and a vector meson, a scalar and a tensor meson, or a pair of vector mesons. 

A main decay file is used (DECAY. DEC) which provides a fairly complete table of the decays 
of particles below the T(4S) to exclusive final states. For generic decays of B mesons, about 50% 
are generated to exclusive states, while the other 50% are produced using Jetset7.4 fragmentation 
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model. Jetset7.4 is also used to generate the cc states and weakly decaying baryons. The original 
Jetset7.4 decay table was modified in order to be in line with recent measurements. 



5.2 The GEANT3 Model 

BBsim is based on the detector description and simulation tool, Geant321 |7j], developed at CERN. 
Geant321 provides tools to construct the detector geometry; to step charged and neutral particles 
through the detector; to simulate the full variety of interactions and decays that each particle species 
may undergo; to register Monte Carlo track hits (referred to as GHits in BABAR terminology); and 
to display the detector components, particle trajectories and track hits. 

BBsim is organized as a set of subsystem packages, each consisting of a standard set of rou- 
tines which are called at various stages of the event simulation. In the initialization phase, the 
detailed subsystem geometries are built from parameters specified in an ASCII geometry database. 
The model includes the definition of the shapes, positions and orientations of all the subdetector 
components, and of the materials used in their construction. 

Long-lived particles ("primaries") produced in the event generation stage are stepped through 
the detector and allowed to interact or decay. In BBsim, hadronic interactions may be simulated 
using either Gheisha @ (default), Fluka @ or GCalor |10fl . Secondaries resulting from interactions 



or decays are tracked as well. However, when recording the Geant3-level genealogy, only the 
kinematics of primaries and secondaries resulting from decays are included. 

Particle GHits are scored in the active detector components. The GHit contains all the infor- 
mation needed to perform the subsequent detector response simulation, such as the position of the 
hit, the direction the track is traveling in, the energy lost, and the Monte Carlo track number. All 
of subsystem GHits are written to an output file, along with the generator- and Geant321-level 
Monte Carlo event generation information. 



5.3 The Detector Response Simulation 

The digitization of track hits is done as a separate step called SimApp. This uses the GHit infor- 
mation as input and outputs digitized data in the same format as the real detector. At the end of 
this stage the Monte Carlo data is processed by the same code that is run on the real data. The 
code is again organized as a set of subsystem packages. These packages contain routines to make 
the simulated data as realistic a match as possible to data taken by the detector. As an example, 
the drift chamber package takes account of the hit efficiencies measured using real data. Several 
figures showing the quality of the agreement between simulation and real data can be found in this 
paper. 

Another function of the code in the SimApp packages is to add in background hits. Rather than 
simulate the backgrounds in the detector, events from samples recorded with random triggers are 
overlayed (using the appropriate luminosity weighting factor) with Monte Carlo simulated events. 

Simulation production takes place at SLAC and at remote sites. The rate of production has 
been increasing with time and currently the level is about three million events a month. Since 
background data is used in the simulation, and the backgrounds in the machine fluctuate, some 
care is needed to ensure that the Monte Carlo is produced with representative backgrounds. Each 
production block corresponds to about a month of real data and uses background files and detector 
efficiencies "typical" of that month. The background events are shuffled to ensure that even small 
Monte Carlo data sets are generated with a sample of background events spread out over the entire 
time period. 
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6 Tracking and Vertexing 



6.1 Track Finding and Kalman Filtering 

Charged track finding in BABAR reconstruction starts with pattern recognition in the DCH. A better 
event time is needed than that available from the L3 trigger, so the first step is to find an improved 
to using three parameter fits (do, cfto, to) to four-hit track segments in the DCH superlayers. 

The first pass at DCH track finding takes the track parameters and associated lists of hits from 
the output of the L3 track finder and performs helix fits; if a fit succeeds, that track is retained. A 
search for additional hits in the DCH that may belong on these tracks is made, and the event time 
is further refined using the timing information of just the hits assigned to tracks. Two subsequent 
track finders are run on the remaining DCH hits. Both track finders work on superlayer segments 
constructed out of hits not used on previous tracks, and are designed to find tracks at lower px 
than found by the L3 trackfinder, to find tracks going through fewer than the full ten superlayers, 
and to find tracks without bias towards the interaction point. In this way more sophisticated track 
finders run in a progressively cleaner tracking environment with a constantly improving measure 
of the event time. 



At the end of this process, the found tracks are refit using a Kalman filter fitter. |11] This fitter 
takes into account the detailed distribution of material in the detector and the full map of the 
non-uniformities in the detector magnetic field. The tracks are then projected into the SVT, and 
silicon strip hits are considered for addition to the tracks if they are consistent with the expected 
error from extrapolating through the intervening material and field. Various possible combinations 
(branches) of SVT hits are considered, and the branch contributing the least to x 2 > while meeting 
pattern requirements designed to limit the number of silicon wafers providing no measurement, is 
retained. When all intersected wafers have been inspected for hits and the best branch is selected, 
a full Kalman fit is performed on the ensemble of DCH and SVT hits. 

Then a search for tracks is performed using the remaining unused silicon strips, as described in 
the next section. At the end of the SVT-only track finding, a search of the track list is made for 
tracks made up of only SVT hits which match closely at the support tube to other tracks made up 
of only DCH hits, since the support tube is a localized source of scattering between the two devices. 
If any candidates are found, a common fit to both sets of hits is attempted. If this succeeds, only 
the combined track is retained; otherwise, both stand-alone tracks are kept on the track list. At 
the end of pattern recognition and fitting, all tracks are checked to see if they lie in momentum 
ranges where the particle hypothesis might effect the result of the fit, and if so, each track is fit 
with up to five stable-particle mass hypotheses. 



6.2 Tracking Efficiency 

The absolute level of track reconstruction efficiency has been determined by studies using control 
samples from data. Changes in the DCH operating point and the reconstruction code during data 
taking are also followed in this way. The polar and azimuthal angle dependence of the track finding 
efficiency is modeled through detailed tuning of the Monte Carlo simulation. The overall level of 
track reconstruction efficiency and its dependence on angle or pt is checked using: 

• Bhabha scattering events at low multiplicity and high pt] events are selected based on 
calorimeter information and the efficiency for finding the e + e~ tracks determined; 

• one-versus-three tau-pair events at moderate multiplicity and px ranging up to about 2 GeV/ c; 
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events are selected without requiring the third track on the 3-prong side, allowing a deter- 
mination of the probability for reconstructing this track with a variety of track selection 
requirements; 

• events containing a D° — > K~ir + ir + ir~ decays where the D° is reconstructed from the decay 
D* + — * D°tt + ; the relative tracking efficiency for tight or loose requirements on the DCH 
portion of the track can be determined, as well as absolute efficiencies using the ratio of 
D° K~ir + TT + n~ to K'tt + ; 

• multiplicity distribution of clean multihadron event selections; 

• comparisons between standalone tracking in the SVT and DCH. 

We assign a conservative 2.5% systematic error per track to the absolute track finding efficiency 
above 200 MeV/ c in px in computing all branching fractions reported at this conference, and treat 
this as fully correlated for all studies that require a DCH portion of the track. This represents the 
level of consistency among track finding efficiency studies done so far, but in no way reflects the 
final systematic error that we expect to achieve with the tracking system. 

6.3 Slow Pion Reconstruction 

The tracking system in BABAR was designed to be highly efficient even for very low momentum 
tracks, which are important for reconstructing decays such as D* + — > D°tt + for which the ir + has 
a very soft momentum spectrum. Tracks with transverse momentum {pr) below about 120MeV/c 
can not be found with the DCH pattern recognition because they curl up and do not penetrate 
deep into the DCH, which has an inner radius is 24 cm. The five-layer SVT was designed to be a 
stand-alone tracking device covering specifically the entire pr range of 50-120 MeV/c. The lower 
Pt limit of 50 MeV/c is the threshold below which pions have a large probability of stopping in the 
material of the beam pipe. 

The SVT stand-alone pattern recognition is performed by two independent algorithms after 
SVT hits have been associated with, and claimed by, stand-alone DCH tracks. The first algorithm 
builds tracks from space points, which are intersections of 4> and z strips. The second algorithm 
forms circle trajectories with the <ft hits and then adds z hits to the circles to form helices. The 
first algorithm is efficient over a wide range of impact parameter and zq values but is susceptible 
to inefficiencies in events with an extremely large number of hits and from silicon modules with 
non- functional sides (4> or z). The second algorithm is less sensitive to extreme combinatorics and 
module defects and is used to find tracks missed by the first algorithm. 

Figure |9| shows the slow pion pr spectra from the decay D* + — > D°tt + , after subtraction of 
combinatorial background using the m(D°) sidebands, for on-resonance data and a mixture of bb 
and cc Monte Carlo which represents the on-resonance conditions. The shape of the px spectra 
agree quite well, which gives us some confidence in our Monte Carlo simulation. Figure ^ also shows 
the tracking efficiency as a function of px as derived from the Monte Carlo. The efficiency turns 
on at a pt of about 50 MeV/c and crosses 80% at about 70 MeV/c. We have tried to measure the 
relative tracking efficiency at low momentum from the data alone using slow tt helicity distributions, 
binned in p*{D* + ). Currently, this measurement is statistically limited, but will eventually yield a 
tracking efficiency measurement at low momentum from the data without Monte Carlo dependence. 
We assign an uncorrelated 2.5% systematic error to the track finding efficiency for those tracks that 
rely on the SVT only. 
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Figure 9: The figure on the left shows the slow pion pr spectrum, after sub- 
traction of combinatorial background using the m(D°) sidebands, for on-resonance 
data (points) and Monte Carlo (boxes). The plot on the right shows the tracking 
efficiency as a function of pr derived from the Monte Carlo. 



6.4 Vertex and Kinematic Fitting 

Vertex and kinematic fitting is widely used to improve four-momenta and position measurements, 
as well as to measure the time difference between decaying B hadrons in the T(4f>) — * BB decay. 
Combined with tracking and particle identification, vertex and kinematic fitting provides the basis 
for event identification. 

The physics properties of a decay are used to apply constraints which translate to better mass 
and position resolutions and larger signal-to-background ratios. For example, in the case of B° — > 
J/ip Kg, the position measurement of the B° can be improved using the fact that the line-of- flight of 
the Kg intersects the J/tp vertex. The energy resolution of the B° can also be improved by applying 
a mass constraint to the J/tp . Several constraints have been implemented utilizing the Lagrange 



Multiplier mechanism [15|: common decay vertex, mass, energy, momentum, beam energy (with 
and without smearing), beam spot position and line-of-flight. In all cases constraints are applied 
in three dimensions. 

One of the fundamental principles in the design was to deal in a simple way with complex decay 
chains. Virtual composite particles and their error matrices are built from the original particles. 
The composite particle then replaces the daughters in subsequent fits and analysis. To account 
for intermediate resonances and neutral particles, vertices can be shared by different virtual states. 
The three-momentum of the virtual particle is fit directly, rather than computed from the updated 
daughters, improving speed and numerical accuracy. 

Non linearities in the fits require the application of an iterative procedure. Simple fits involving 
only vertex constraints (except long-lived particles, called V°'s) are, however, accurate enough with 
a single iteration. This has an important impact on the amount of time consumed in vertex fitting. 
The other fits involving kinematic constraints and V 0, s require, in general, more than one iteration. 
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The criteria we use for convergence is that the change in x 2 between two successive iterations is 
less than 1%, with a maximum of six iterations. In all cases the initial location of the vertex is 
estimated by solving analytically (with an iterative procedure) the point of closest approach of the 
two tracks, using a second-order approximation at each point. 

Track and vertex parameter errors depend heavily on our understanding internal and global 
alignment, detector material, and magnetic field variation. A set of control samples (77 — > 4tt, 
t — ► 3-prongs, 'Tr + ir 7r + , B — > J/ip K + ) has been used to estimate from data the 

vertex position and mass resolutions and systematic uncertainties in their determination. The pull 
distributions of these control samples show that tracking and vertex errors are estimated correctly 
at the 10-20% level and that biases in determining vertex positions due to alignment problems are 
at the level of a few microns. 

6.4.1 Beam spot and primary vertex determination 

In order to accommodate the movements of the PEP-II beam with respect to the BABAR detector, 
the position of the collision point at PEP-II is determined from track-based methods. The apparent 
size from two-prong events is about 150 /jm in x, 50 //m in y and 1 cm in z. The apparent size in the 
y direction is totally dominated by the track resolution. A better estimate of a y can be obtained 
from the knowledge of the luminosity, the beam currents and the size in x, providing a value of 
about 4 /j,m, varying within 10% on a time scale of hours. 

For many physics analyses and reconstruction tasks, it is necessary to have a more accurate 
estimate of the collision point in x and z. For instance, reconstruction of 7 and ir° candidates uses, 
in many cases, the primary vertex as the point of origin in reconstructing vectors. The primary 
vertex^ is estimated on an event-by-event basis from a vertex fit which uses charged tracks with 
an impact parameter (with respect to the beam spot position) less than 1 mm in the transverse 
plane. Tracks with high x 2 contribution to the vertex fit are removed until an overall x 2 probability 
greater than 1% is obtained. The resolution is about 70 //m in x, y and z for hadronic events. 

6.4.2 Reconstruction of the tagging-!? vertex 

In measurements of the B lifetime, mixing rate and time-dependent CP-violating asymmetries, 
we fully reconstruct or partially reconstruct only one B and then determine the distance between 
the two B decays. Determination of the fully-reconstructed B decay vertex makes use of all the 
final-state tracks of the B decay chain. The other B vertex (tagging- i3 vertex) is reconstructed 
with the remaining charged tracks in the event. To retain high efficiency this is done with inclusive 
techniques. Complications arise from the fact that one has to deal with secondary tracks from 
short- and long-lived particles. To deal with these problems and minimize their impact, the three- 
momentum of the tagging- B (called a pseudo-track) and associated error matrix are derived from 
the three-momentum, decay vertex and error matrix of the fully reconstructed B candidate, and 
from a knowledge of the average position of the interaction point and the T(4S) four-momentum. In 
the case of partial B reconstruction (where the three-momentum is not known) the constraint from 
the B pseudo-track cannot be used. This pseudo-track is fit to a common vertex with the remaining 
tracks in the event. Reconstructed V° candidates are also used in this procedure, reducing biases 
from long-lived decay particles. Tracks with contributions to the x 2 greater than six are removed 

For light quark events, the reconstructed vertex is a good estimate of the true primary vertex. For BB events, 
where both B mesons travel along the z axis in the laboratory frame, the primary vertex provides an average B decay 
position. 
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from the fit. The procedure is iterated until there are no tracks contributing more than six units 
to the x 2 or an tracks are removed. 

The Az resolution function obtained with this technique can be accurately described by the 
sum of two Gaussian distributions with different means and widths. Typically, 80% of the events 
are in the core Gaussian, which has a width ~ 100 [im. The remaining events reside in the tail 
Gaussian which has a typical width ~ 350 /xm. The bias due to secondary tracks is ~ — 20/iin for 
the core Gaussian and ~ —80 /xm for the tail Gaussian. An inclusive D*~ — > D°ir~ control sample 
(selected from continuum and off-resonance running) is used to check the Az reconstruction and 
resolution for data. 



7 Particle Identification 

7.1 Charged Particle Identification 
7.1.1 Electrons 

Electrons are identified with the EMC, the DCH and the DIRC systems, using quantities such as the 
ratio of the deposited energy over the momentum, and variables describing the spatial shower de- 



velopment. The lateral shower development is measured with a method introduced by ARGUS [16| 
and expanded into various moments containing information about the shower's azimuthal distribu- 
tion with respect to the particle's initial momentum [|7]]. A cut is applied on the energy loss using 
a truncated mean of 40 samples (maximum) recorded in the DCH. The measured Cherenkov angle 
in the DIRC is required to be within 3a of the expected value for the electron hypothesis. The 
efficiency of electron identification is studied by using electrons obtained from radiative Bhabhas 
and events of the process 77 — > e + e~ . The misidentification probability for pions is measured with 
3-prong t decays. The tight selector has an average efficiency of 94.8% in the momentum range 
0.5 GeV/c < p < 2 GeV/c, with a pion misidentification probability of roughly 1.2%. The very tight 
selector has an efficiency of 88.1% with an average pion misidentification of 0.3%. The dependence 



of the efficiency on momentum and polar angle is shown in Figure 10 



7.1.2 Muons 

Charged tracks reconstructed in the DCH are extrapolated to the flux return iron using a detailed 
map of the non-uniform magnetic field and accounting for the expected average energy loss. The 
predicted average position of the intersections with the active detector planes is computed, including 
the uncertainty due to multiple scattering. All hits found in each readout view within a specified 
maximum distance from the predicted intersection are associated with the charged track. We 
require, for tracks within the acceptance of the EMC, the energy deposition to be consistent with 
that expected of a minumum ionizing particle. A signal in at least two layers of the IFR is also 
required. A cut is also made on the difference between measured and predicted (based on the 
muon hypothesis) total number of interaction lengths traversed in all subdetectors. We expect the 
average number of signal strips per layer to be larger for pions produced in an hadronic interaction 
than for muons. The average value and the r.m.s. of the strip pattern for the different hit layers 
provides fi/ir discriminating power. Finally, we reject tracks with a large track-match x 2 or a large 
X 2 for a polynomial fit to the IFR cluster. 

The performance of the muon selectors has been tested using control samples obtained from 
data. The processes considered are /U/xee and /z/ry for muons, and r — > 3-prongs and K$ — ► 
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Figure 10: Electron identification efficiency and pion misidentification for the very tight selection 
criteria as a function of momentum (left) and polar angle (right) (note the different scales for 
electrons and pions). 
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Figure 11: Muon identification efficiencies (left) and pion misidentification rates (right) for the 
tight muon selector, showing average efficiencies of 75% and misidentification rates of 2.5% for the 
momentum range 1.5 to 3GeV/c. 

7r + 7r~ for pions. The selection of the control samples is based on kinematic variables and particle 
identification. No muon requirements are applied to the tracks that are used as a control sample, 
in order to avoid any bias in the determination of the particle identification performance. Different 
muon selection criteria with different levels of purity are defined for analysis specific applications. 
In the range 1 < p < 3GeV/c the average efficiency is about 77% for muons with a 2.5% pion fake 
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rate (see Figure 



7.1.3 Charged hadrons 

Kaons are selected using the information from the SVT, DCH, and DIRC systems. Likelihoods are 
calculated from the dE/dx measurements in the SVT and DCH, and from the angle and number 
of photons found in the Cherenkov ring in the DIRC. The individual likelihoods from the three 
detectors are then multiplied together and a cut on the ratio of likelihoods (K vs tt, K vs p) is 
applied: 

L n • r < Lk and L p • r < Lk (4) 

Various kaon selections have been implemented, achieving different levels of efficiency and purity. 
Each detector output is used only in a specific interval in momentum. For example, for a tight 
selection, the ranges used for each subdetector are 0.025-0.7 GeV/c for SVT, 0.090-0.7 GeV/c for 
DCH, and > 0.6GeV/c for the DIRC. The minimum requirements for a track to be considered 
in the individual subdetectors are 2 and 10 dE/dx samples for the SVT and DCH respectively, 
and track extrapolation within the DIRC acceptance (number of expected photons greater than 
zero). Finally, the value of likelihood ratio r was dependent on the momentum of the tracks to be 
identified. In the tight mode, the intervals and ratio requirements were: 

0.5-0.7GeV/c r = 15 
0.7-2.7GeV/c r = l 
>2.7GeV/c r = 80 

The performance of the various selections has been measured with data using control samples 
of kaons and pions from the kinematically identified decays: K s — > -k + it~ and D* + — > D°(^ 
K~Tr + )n + . To select the K® sample, very tight cuts are used (e.g. on the angle between the 
Kg candidate vertex direction and the vertex momentum direction, on the decay distance from 
the interaction point, on the reconstructed Kg candidate mass). The plot of the invariant mass 
of the tt + tt~ pairs is shown in Figure [T^(left); the purity of the resulting pion sample is > 99%. 
The D* + sample, obtained with a tight cut on the difference between the Ktttt and Kn masses 
(0.1445 < AM < 0.1465 MeV/c 2 ) shown in Figure [i~2](right), has a combinatorial background of 
around ~ 13% (for a kaon sample purity of ~ 90%, with the purity higher at higher momenta). In 
determining particle identification performance, a correction for the background contribution has 
been made in each momentum bin using sideband subtraction. 



The measured kaon efficiency and pion misidentification rates are shown in Figure 13 for the 
tight selection based the D* + sample. The efficiency is found to be 90%, while misidentification 
rate is 2.5%, averaged over the full momentum range. Note that these results are highly dependent 
on the momentum spectrum of the kaons and pions in the particular samples used in an analysis. 



7.2 Neutral Particle Identification 
7.2.1 Photons 

Photons are detected in the EMC where adjacent crystals with energy deposits exceeding 1 MeV 
are grouped in clusters. Channels marked as noisy by online or offline monitoring are excluded. 
Clusters with more than one local energy maximum are then split into "bumps" and the energy 
of each crystal is partially assigned to each bump by a simultaneous iterative adjustment of the 
centers and energies of the bumps, asuming electromagnetic shower shapes. In the next stage all 
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Figure 12: Mass distribution (left) of Kg candidates used to select a pion control sample. Distri- 
bution (right) of the difference of the D* and D° candidate masses, used to select control samples 
of kaons and pions. 
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Figure 13: Kaon efficiency (filled points) and pion misidentification (open points) for a tight kaon 
selection, measured with a D* decay control sample. 

charged tracks reconstructed in the tracking volume are extrapolated to the EMC entrance and 
a track-bump matching probability is calculated for each pair. All bumps with a small matching 
probability are treated as photon candidates. A small number of these candidates, where the 
bump shape is incompatible with that expected of an electromagnetic shower, are rejected. In 
the majority of physics analyses, photon candidates lying anywhere within the instrumented EMC 
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volume are used if their energy is above 30MeV. 



7.2.2 7T° and rj Reconstruction 

Neutral pion and r/ candidates are formed from pairs of photon candidates assumed to originate 
from the interaction point. The invariant mass spectrum of all such pairs is shown in Figure [l4| for 
(E 7 >30MeV, E 77 >300MeV) and (E 7 >100MeV, E 77 >lGeV), where clear vr° and r/ peaks can 
be seen. The tt° mass resolution in multihadron events is 6.9MeV/c 2 , while in low occupancy tt 
events it is 6.5MeV/c 2 for tt° energies below about 1 GeV/c. Improved resolution is also observed in 
hadronic events where only isolated photons are considered. The r\ mass resolution is 16MeV/c 2 . 



The tt and i] mass peaks from inclusive multihadron events are also shown in Figure |T4[ The signal- 
to-background ratio for the r\ can be improved by vetoing photons participating in candidates in 
the 7T° signal region. Combinatorial background in hadronic events is supressed by a factor of 5 by 
this veto requirement, with a signal loss of only 25%. 




3500 
3000 - 
2500 
2000 - 
1500 
1000 
500 - 



0.1 0.2 0.3 0.4 

77 mass 



0.5 0.6 0.7 0.J 

(GeV/c 2 ) 




0.4 0.45 0.5 0.55 0.6 0.65 0.7 

77 mass (GeV/ c 2 ) 



Figure 14: Invariant mass of 77 pairs from hadronic events. Upper pair of plots are for E 7 >30 MeV 
and E 77 >300MeV, with enlargement of the vt° mass region on the right; lower pair of plots are 
for E 7 > 100 MeV and E 77 >1 GeV, with enlargement of the 77 mass region on the right. 

The detector segmentation and achieved spatial resolution allow the reconstruction of tt°'s 
with photon separation in the electromagnetic calorimeter as small as 5 cm, although with some 
deterioration of the resolution. The small fraction of high energy ir° whose photons cannot be 
resolved into individual photon bumps (~10% in the 4-6 GeV region) can be separated from single 
photons using the cluster shape (the second moment of the energy distribution around the cluster 
centroid) . 

Large samples of inclusive and exclusive data have been used for simulation comparison studies. 
The 7T° mass value, resolution, and efficiency, and their (weak) energy dependence were found to 
be in very good agreement. Furthermore, the ir° mass, width, and yield, are monitored in real time 
during data collection and reconstruction. Limits on the possible discrepancies between the data 
and the simulation for the photon energy scale, resolution, and reconstruction efficiency, have been 
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conservatively set to be 0.75%, 1.5% and 2.5% respectively. 



7.2.3 Neutral Hadrons 



K° s -» vr+vr- 

To reconstruct K s candidates, we consider any pair of oppositely-charged tracks combined in 
a common vertex. No cut is made on the probability of x 2 ■> but the vertex fit must converge. 
The recontructed momentum of the Kg candidate is required to be aligned with the direction 
between the interaction point and the reconstructed vertex. Finally, a minimum transverse 
momentum of the daughters with respect to the flight direction is required in order to eliminate 



combinatorial and A contamination. For Figure 15, a tighter selection is applied on a small 
data sample; the x 2 probability is required to be greater than 0.01, the angle between the 
direction of flight and the direction from the interaction point to the reconstructed vertex is 
required to be greater than 45 mr, the transverse momentum of the daughters with respect 
to the Kg flight direction is required to be between 0.11 and 0.22 GeV/c. 
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Figure 15: K s mass resolution (tight selection). The fit is performed with a double Gaussian, the 
central Gaussian contains 70% of the events and has a width of 2.82 ± 0.19MeV/c 2 . 



K° s 7r°7r° 



Non-overlapping 7r° candidates are combined to construct K® — > candidates. For each 
candidate with an energy above 800 MeV and a mass between 300 and 700MeV/c 2 , we de- 
termine the most probable decay vertex along the path defined by its momentum vector and 
the primary vertex. The point where the product of the probabilities from the two ir° mass 
constraint fits is a maximum is chosen as the K® decay vertex. We demand that this point 
lie in a region between —10 and +40 cm from the primary vertex and that the K® mass at 
that point is in the range 470-536 MeV/ c 2 . This method |§ improves significantly the mass 
scale and resolution. A vertex resolution of 5 cm along the flight direction is predicted by the 
simulation. The inclusive K® mass peak from real data, exhibiting a resolution better than 
10 MeV/ c 2 , is shown in Figure 16. 
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Figure 16: ir ir mass spectrum and K® signal from hadronic events 

The EMC and IFR detectors are both used for the detection of Kl mesons, which are recon- 
structed as neutral clusters that cannot be associated with any charged track in the event. In 
the EMC, the background consists primarily of photons and high energy 7r°s. We are able to 
distinguish Kl clusters from these backgrounds by using a shower shape analysis, performed 
by expanding the shape into a series of Zernike moments. For clusters with an energy greater 
than 500 MeV, a single Zernike moment (2,0) is sufficient to discriminate Kl mesons and 
photons. Between 200 and 500 MeV, a two-dimensional cut ((2,0) and (4,2)) is used. Below 
200 MeV, we are unable to sufficiently discriminate between photons and Kl mesons. IFR 
clusters are defined and selected as hits in two or more resistive plate chamber layers. The 
background is dominantly charged particles and detector noise. Some splitoffs from charged 
hadronic showers are missed by the tracking association, due to the irregular structure of 
these showers. We suppress these clusters by rejecting Kl candidate clusters close to any 
track in the event. 

The performance of the Kl identification is demonstrated in Figure |l7|a, using events from 
e + e~ — > 7(0 — > KsKl)- In these events, the position and energy of the Kl can be constrained 
from the measured photon and K$ kinematics. Figures |l~7| b and |i~7| c illustrate the performance 
of the Kl selection. We observe good agreement between Monte Carlo simulations and the 
data. The efficiency for Kl reconstruction agrees to within the current statistical error, which 
is approximately 10% for both IFR and EMC selections. 
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Figure 17: The missing mass distribution (p 2 K ) of candidate <fyy events is shown in (a). We find 
594 ± 40 e + e~ — > j((f) — > K$Kl) events. Using these events, (6) shows the comparison of the 
Zernike moment (2,0) distribution for clusters from Ki mesons (points) and photons (histogram). 
A comparison between data (points) and Monte Carlo (histogram) is shown in (c) for the number 
of hit layers in the Instrumented Flux Detector (plot ( c) is absolutely normalized) . 



8 Common Analysis Issues 

8.1 Centre of Mass Energy and Boost 

The beam energies provided by PEP-II are obtained with an algorithm capable of tracking variations 
to within about 1 MeV, and are used to determine the time-dependence of the centre of mass energy. 
The absolute energy scale is determined from a sample of fully-reconstructed B candidates (a recent 
measurement of the B mass by CLEO [13] is used). 
We define two main coordinate systems: 

• A BABAR coordinate system (LAB), linked to the BABAR detector. The z axis is along the 
DCH axial wires, the y axis is vertical (upward) and the x axis points towards the exterior of 
the PEP-II rings. The BABAR z axis and the PEP-II electron beam are not exactly aligned. 
The former is tilted in the horizontal plane by 19mrad with respect to the beam axis. 

• A centre-of-mass coordinate system (CMS), a frame where the two beams have equal energy 

= \/s/2 (where is the centre-of-mass beam energy). The z axis of the CMS lies 
along the electron beam direction (this is the relevant choice for most physics quantities). 
The Lorentz transformation from LAB into CMS quantities is the product of a rotation that 
aligns the z' axis of the rotated frame with the boost direction, and a Lorentz boost along the 
new z' into the CMS frame. The rotation angles are determined on a run- by-run basis from 
the opening angle and flight direction of fi + fi~ and e + e~ events, and are verified against the 
orientation of the collision spot for the beams. The magnitude of the boost is obtained from 
the PEP-II energies. 



8.2 Continuum Rejection Variables 

We use common shape variables, computed with charged tracks in the centre of mass frame, to 
preferenctially reject continuum events. The most important of these variables are thrust, sphericity 
and R2, the ratio of 2 nd to th Fox- Wolfram moments. 
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• Thrust: The thrust axis, T, of an event is defined to be the direction that maximizes the 
sum of the longitudinal momenta of the particles. The value of the thrust T must lie in the 
range 0.5 to 1, where T ~ 1 corresponds to a highly directional event and T ~ 0.5 to an 
isotropic event. 

• Sphericity: Sphericity is a measure of the sum of squares of transverse momenta for each 
track with respect to the event axis. Highly directional events have low sphericity whereas 
isotropic events tend to have sphericities close to one. 

• Fox- Wolfram moments: The £ til Fox- Wolfram moment is momentum-weighted sum of 
Legendre polyomial of the £ th order computed from the cosine of the angle between all pairs 
of tracks. The ratio R2 of Fox- Wolfram 2 nd to th moments is the variable which provides the 
best separation between BB signal and continuum. Jet-like continuum events tend to have 
higher values of R2 than the more spherical BB events. 

8.3 Event Selection and B Counting 

8.3.1 Fiducial regions 

The acceptance region in polar angle for charged tracks, corresponding to full SVT coverage, ranges 
from a minimum of 410 mr with respect to the z axis in the forward direction to 602 mr in the 
backward direction (0.41 < 9 < 2.54). For neutral clusters, the fiducial region is reduced to an 
angle 732 mr with respect to the z axis in the backward direction (0.41 < 9 < 2.41) to account for 
the more limited EMC coverage in the backward region. 

8.3.2 Multi-Hadron selection 

Multi-hadron event selection is designed to have a very high efficiency for BB events while keeping 
the systematic uncertainties on the determination of that efficiency as low as possible. The contri- 
bution of continuum and r pair events in the sample is reduced by a cut on R2. Using only tracks 
and clusters in fiducial volume described above, the selection criteria are as follows: 

• Four or more charged tracks, at least three of which must have associated DCH information 

• Primary vertex within 5 mm in x and y of the nominal beam spot position, 

• R2 less than 0.7, 

• Total energy greater than 5 GeV. 

The efficiency for BB events is 96%, while the contamination due to beam gas, two-photon and 
tau pair events is approximately 2%. 

8.3.3 Muon pair selection 

In the centre of mass frame, events containing two charged tracks with momenta larger than 2 GeV/c 
and 4GeV/c, respectively, and with angle with respect to the beam axis greater than 725 mr are 
selected. The two tracks must have an acolinearity smaller than 10°. The mass of the pair must be 
larger than 7.5 GeV/ c 2 . The total energy deposited in the calorimeter by the two particles must be 
smaller than 1 GeV. The event is rejected if neither track has an associated energy measurement. 
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8.4 Offline Luminosity Determination 



The integrated luminosity recorded by BABAR is determined offline using Bhabha, muon pair, and 
gamma-gamma events. The measurement precision is limited by systematic errors in all cases to 
the several percent level, although this clearly can be improved to 0.5% or better. The systematic 
uncertainty for the Bhabha measurement is dominated by the efficiency corrections for the Bhabha 
veto by the Level 3 trigger, and by theoretical uncertainties on the differential Bhabha cross section. 
Likewise, we are still in the process of comparing and understanding Monte Carlo generators for 
muon pair events, including proper treatment of radiative events. A better understanding of the 
impact of occasional hot towers in the EMC will improve the systematic error on the gamma-gamma 
luminosity. The luminosity determinations are consistent at the several percent level, but are much 
better than this for tracking relative luminosity and the variation with beam conditions. 



8.5 B Counting 

We estimate the number of BB pairs N BB - from the total number of events passing the hadronic 
event selection Nmh and the total number of events passing the muon pair selection N^: 

N° n 

N bS = N^ h -kxN^ h x-^ , (5) 

where the superscript on refers to on-resonance data and off to off-resonance, k is a factor close to 
1 that corrects for changes in efficiency and cross section with center of mass energy. It is estimated 
by Monte Carlo to be 0.9962. The systematic error on N BB is 1.7%, dominated by the variation 
between running periods of the ratio Nffi H /N°l/ _. 

Branching fraction measurements frequently use the number of produced BB pairs. The effi- 
ciency for BB events to pass the selection criteria is determined by Monte Carlo to be 96.0%. The 
total systematic error on the number of produced BB pairs is 3.6%. 



8.6 Exclusive B Reconstruction Variables 
8.6.1 Kinematic variables 

A number of kinematic quantities could be used to characterize the reaction e + e~ — > BB where one 
of the B mesons is fully reconstructed. However, our goal is to select a pair of such kinematic vari- 
ables, having little correlation between them, that make use of the maximum available information 
for optimal background rejection. 

At the T(4S) resonance the B mesons are produced with very small Q values. As a result their 
center-of-mass energy E* and momentum p* are very sensitive to fluctuations of the center-of-mass 
beam energy E£ = y/s/2. In contrast, the variable AE defined as: 

AE = E* — El , (6) 

is relatively insensitive to E% fluctuations. 

The distribution of AE is peaked at zero for e + e~ — ► BB events, and its width is governed in 
most cases by the beam energy measurement error a\. It is not necessary to boost the B candidate 
to the center-of-mass frame to compute AE, as can be seen by writing AE in an explicitely Lorentz 
invariant form: 
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where P and Pj are the Lorentz vectors representing the B candidate four momentum and the 
initial-state four momentum respectively. 

The AE variable is used in conjunction with either one of the following two mass variables: 



• Beam-energy substituted mass mES) defined as: 



mES 



±s + gVgj 



P 2 ■ (8) 



where E^ and are the total energy and the three momentum of the initial state in the 
laboratory frame, and is the three momentum of the B candidate in the same frame. At 
a symmetric machine (where Jf i = 0) this variable is obtained simply by substituting the 
beam energy E£ from the measured B candidate energy E*: 



™ES = V E f ~ P* 2 ■ ( 9 ) 

This is the definition of the variable used in ARGUS and CLEO publications under the 
name "beam-constrained" mass. The advantage of using definition of Eq. || at an asymmetric 
machine, where one needs to assign masses to candidates in order to boost to the center of 
mass frame, is that tties is computed in the laboratory frame without any prior commitment 
to the identification of particles among the B daughters. 

Beam-energy constrained mass mEC> defined as: 



mEC = V^ 2 -^ 2 » ( 10 ) 

where E and p are center-of-mass frame quantities obtained by performing a kinematic fit 
with the E* = E^ constraint. 

The two choices (AE, toes) an d (AE, toec) are almost equivalent. For certain analyses, it may 
be helpful to exploit the property that ?ties does not require daughter mass assignments, while the 
kinematical fit performed for toec makes optimal use of the information collected by the detector, 
in addition, the pair (AE, rriEc) has the smallest correlation. 

8.6.2 Background parametrisation 

To describe the background shape in beam-energy constrained mass plots, we use the ARGUS 
distribution [jlj] given by: 

6(m < mo) / ' ; ~Z r 2 1 

A(m\mo, c) = — x my 1 — (m/mo) x exp c(l — m/mo) , (11) 

where tuq represents the kinematic upper limit and is generally held fixed at the center-of-mass 
beam energy E£. 



43 



9 Conclusions and Prospects 

Since the start of running at PEP-II in May 1999, the BABAR experiment has accumulated over 
12.7 fb -1 of data. Many of the detector sub-systems have reached their design goals and significant 
progress has been made towards understanding the performance of the detector. Several preliminary 
measurements with the BABAR detector are to be presented at this conference. We anticipate having 
a data set equivalent to a significantly larger integrated luminosity by the end of this year. Thus, 
the results presented here are just the beginning of an exciting program of physics with BABAR at 
PEP-II. In future years, we expect to be able to undertake detailed studies of CP violation in the 
B meson system. 
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